Hs27 fibroblast response to contact guidance cues

Contact guidance is the phenomena of how cells respond to the topography of their external environment. The morphological and dynamic cell responses are strongly influenced by topographic features such as lateral and vertical dimensions, namely, ridge and groove widths and groove depth (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\text{R}}}_{{\text{w}}}, {{\text{G}}}_{{\text{w}}},\,\mathrm{and}\,\,{{\text{G}}}_{{\text{D}}}$$\end{document}Rw,Gw,andGD, respectively). However, experimental studies that independently quantify the effect of the individual dimensions as well as their coupling on cellular function are still limited. In this work, we perform extensive parametric studies in the dimensional space—well beyond the previously studied range in the literature—to explore topographical effects on morphology and migration of Hs27 fibroblasts via static and dynamic analyses of live cell images. Our static analysis reveals that the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\text{G}}}_{{\text{D}}}$$\end{document}GD is most significant, followed by the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\text{R}}}_{{\text{w}}}$$\end{document}Rw. The fibroblasts appear to be more elongated and aligned in the groove direction as the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\text{G}}}_{{\text{D}}}$$\end{document}GD increases, but their trend changes after 725 nm. Interestingly, the cell shape and alignment show a very strong correlation regardless of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\text{G}}}_{{\text{D}}}$$\end{document}GD. Our dynamic analysis confirms that directional cell migration is also strongly influenced by the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\text{G}}}_{{\text{D}}}$$\end{document}GD, while the effect of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\text{R}}}_{{\text{w}}}$$\end{document}Rw and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\text{G}}}_{{\text{w}}}$$\end{document}Gw is statistically insignificant. Directional cell migration, as observed in the static cell behavior, shows the statistically significant transition when the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\text{G}}}_{{\text{D}}}$$\end{document}GD is 725 nm, showing the intimate links between cell morphology and migration. We propose possible scenarios to offer mechanistic explanations of the observed cell behavior.

the intrinsic morphology and/or variation of respective cell type in the intracellular population are closely related with their physicochemical functions 26 .Many studies have simultaneously varied two dimensions, R w and G w , while poorly controlling over the other ( G D ), which can lead to inaccuracy in predicting cell behaviors because of inconclusive data generated by intertwined topographic effects.
Despite the inconsistency, it has been well accepted that contact guidance technologies considerably alter cell morphology 11,27,28 and migration 25,28 and, as a result, have been considered as a possible mechanism to enhance wound healing processes, e.g., by guiding cells to wound sites and controlling their elongation, alignment, and directionality.In this regard, fibroblasts become relevant due to their important roles in wound healing process including breaking down fibrin clots and repairing ECM and collagen structures 29 .However, experimental characterizations of anisotropic morphological behaviors of fibroblasts on contact guidance surfaces are still limited to a relatively small range of topographical dimensions, i.e., G w and R w ≤ 2 µm or G D ≤ 300 nm 11,25,28 .Furthermore, the effect of contact guidance on dynamic cell behavior has been rarely performed and, to our best knowledge, the available results are limited to G w and R w ≤ 2 µm with G D ≤ 300 nm and G w and R w = 3, 5, and 10 µm with G D ≤ 300 nm 25,28 .However, considering typical dimensions (diameter > 10 µm and thickness ∼ 550 nm ) of fibroblasts 30 , it is worth investigating cell anisotropic behaviors on the substrate having intermediate feature sizes of widths and depths ( 2 µm ≤ G w and R w , and 330 nm ≤ G D ≤ 1000 nm ) since each fibroblast can spread on enough number of topographical cues to impose constraints along the groove direction, which is sufficient to affect their morphology and migration.
Beyond the necessity of a contact guidance study for fibroblasts with a broad range of topographical dimensions, it would be highly desirable to independently evaluate a contribution for individual topographic features.To our best knowledge, no effort has yet been made to perform a parametric study for fibroblasts in order to reveal each topographical effect, while holding one-dimension constant with varying the others.The quantitative characterization of cellular behaviors under influence of well controlled yet diverse combinations of geometrical cues is still challenging due to lack of a fabrication technique that allows precise control of small-scale topographic features and multiplexed cell studies under identical extracellular environments.Without such unique capabilities, characterizing the correlation between individual effects of topographical cues and the corresponding cellular behavior with statistical significance becomes non-trivial.On top of that, several hypotheses have been proposed to explain their underling mechanisms; some studies suggest that mechanical interlocking or capillary force results from the groove of the substrates 31,32 and others present that cells lying on the micro-or nanotextured surfaces are able to reorganize their fibers network for balancing internal and external forces 33,34 .Characterizing the specific effects of individual topographical features, using a time-effective high throughput method, is crucial to explore possible mechanisms for cell-type dependent anisotropic cell behaviors and to provide effective design space for biomedical applications of contact guidance.
In here, we characterize the morphological and dynamic responses of Hs27 fibroblasts on multiplexed monolithic quartz chips 35 with a wide range of topographical dimensions well beyond what has been previously reported.Our multiplexed approach allows quantitative characterizations of R w , G w , and G D dependent cell responses while minimizing unwanted biological variations during cell culture.First, we perform static analyses on three different morphological measurements, such as cell spread area, aspect ratio and cell alignment degree, due to their biological and biomedical relevance.Then, we quantify directional cell migration in the scope of effectively guided cell motion.Finally, we utilize both the cell membrane deformation and focal adhesion models to offer mechanistic explanations for the observed cell behavior of HS27 fibroblasts on the multiplexed contact guidance chips.

Fabrication of platforms
To examine the effects of substrate anisotropy on cellular alignment, we created multiplexed patterns with an emphasis on achieving uniform smoothness on the entire substrate including the patterns.The surface topographical cues consist of ridge and groove widths, and groove depths using deep etching process and seeded Hs27 fibroblasts on the etched substrates as shown in Fig. 1a.To be specific, a 20 nm Cr thin film was e-beam evaporated onto 25 mm diameter fused silica coverslips and subsequently coated with AZ1518 photoresist (Microchemicals GmbH) and baked on a vacuum hot plate at 100 °C for 90 s.The photoresist was patterned in a Heidelberg VPG200++ laser pattern generator at a dose of 60 µC/cm 2 and developed with the use of using AZ developer (Microchemicals GmbH) and rinsed in DI water.The patterned Cr film acted as hard mask, allowing only the exposed fused silica area to be etched (Oxford PlasmaLab 100 ICP380 system) by CHF 3 and Ar, with flow rates of 10 and 15 sccm, respectively, at 3 mTorr and 20 °C to achieve various topographical conditions (more details are given in Ref. 35 ).To systematically control topographical cues, ridge and groove widths ( R w and G w , respectively) ranging from 2 to 10 μm with different R w -to-G w ratios were cofabricated on each multiplexed contact guidance chip.Three different groove depths (i.e., G D = 330, 725, and 1000 nm) were used to quantify the effect of groove depth on cell response.The detailed dimensions of the contact guidance chips are presented in Table S1.The lower limit of R w = 2 µm was motivated by the size of typical mature FAs 36,37 , whereas the upper limit of ridge or groove width is 10 μm, which is large enough not to confine focal adhesions and is also about the diameter of the nucleus 38 .

Cell culture
For cell culture, Hs27 cells (ATCC, CRL-1634, Lot #7006836) were cultured in DMEM (ATCC, 30-2002) with 10% FBS (ATCC, 30-2020, Lot #802500) at 37 °C and 5% CO 2 as per manufacturer instructions, without antibiotics.The contact guidance chips and surrounding chambers were assembled under aseptic conditions.Ethanol was introduced to the chambers for 5 min and then N 2 dried for additional sterilization of the contact guidance chips.FN at 25 µg/mL in 10 mM PBS was drop coated directly onto the chips and allowed to adsorb for 1 h at room temperature, and subsequently rinsed 5× with DMEM to remove any excess FN.The contact guidance chambers were then allowed to equilibrate Zeiss microscopes featuring environmental chambers at 37 °C, 5% CO 2 and 95% humidity for at least 30 min before the introduction of cells.Hs27 cells were harvested for contact guidance experiments between 30 and 60% confluence in their logarithmic growth phase.Complete media was aspirated and replaced with 0.25% Trypsin-EDTA (Gibco, #25200-056) for 5 min at 37 °C, after which serum free media was introduced.

Data analysis
The morphology of Hs27 fibroblasts cultured on the monolithic quartz platform was evaluated by analyzing the shapes of all cells adhering to patterned geographies.Image analysis was conducted under 10× phase microscopy and those images were acquired every 10 min and converted into TIFF files using Zen software.For example, the optical images, taken 24 h after initial cell seeding using a 10× phase objective, in Fig. 1b shows Hs27 cells on G10R10 ( G w = R w = 10 µm ) with G D = 330, 725, and 1000 nm from left to right, respectively.The images show that the cells exhibit clearly different cell shapes depending on G D .By using ImageJ (version 1.41, NIH), we manually draw the outline of each cell on the substrates) and then, quantitatively measured the area of individual cells (see the Supplementary Fig. S1).The cell outline was then fitted using an ellipse to estimate the cell aspect ratio ( α ) using a length ratio of the major to minor axes and then obtain the cell alignment angle ( θ ) between the major axis and the direction of groove patterns (see Fig. 1a).
The cell's time-dependent coordinates were measured by tracking its nucleus via ManualTracking plugin in the ImageJ between 10 min time frames.For dynamic analysis, we only consider the following case: cells migrate .The outline of a cell is fitted by an ellipse to estimate a length ratio of the major to minor axes and a cell alignment angle ( θ ) between the major axis and the direction of the groove patterns.(b) Optical images of live Hs27 cells on G w = R w = 10 µ m (i.e., G10R10) with (left) G D = 330 nm, (middle) = 725 nm, and (right) = 1000 nm, respectively.more than 6 h and their nucleus are on the topographical structures, followed by Ref. 35 .This data has been used to compute several dynamic measurements (directional orientation, angular displacement, speed, mean square of distance, directionality ratio, and velocity autocorrelation).For all image analyses, we characterized more than 20 cells for each of the topographical dimensions including control experiments on a flat substrate.For precise analysis of cellular static behaviors, we ruled out the following two cases: (1) cells underwent mitosis and (2) cells spread over both flat and patterned areas of a substrate.Statistical analysis was performed by Kruskal Wallis nonparametric test and Dunn-Sidak post hoc test using MATLAB for the estimation of significance level and statistical comparison of each pair of groups.All data figures were prepared by MATLAB where the central mark indicates the median, and the bottom and top edges of the box point to the 25% and 75%, respectively.The outliers are plotted individually if they are more than q 3 + 1.5 • q 3 − q 1 or less than q 1 − 1.5 • q 3 − q 1 , indicated by the '+' symbol in which q 1 and q 3 are the 25th and 75th percentiles of the data set, respectively.

Results
We experimentally characterize cell morphology and migration of Hs27 fibroblast on the monolithic quartz contact guidance platform.Our goal is to systematically quantify the effect of lateral (ridge width ( R w ), groove width ( G w )), and vertical (groove depth ( G D )) dimensions on static and dynamic behavior of fibroblast cells.For the lateral dimensions, we consider three different cases, namely, G2NR, GNR2, and GNRN.First, G2RN indicates fixed groove width (i.e., G w = 2 µm ) with different combinations of ridge width (i.e., R w = 2, 3, 4, 6, 8, and 10 µm ).Similarly, GNR2 corresponds to R w = 2 µm with different G w values (i.e., G w = 2, 3, 4, 6, and 8 µm ).G2RN and GNR2 are designed for independent characterization of the R w -and G w -dependent cell responses, respectively.Finally, GNRN represents topographical cues with R w /G w = 1 where R w = 2, 3, 4, 6, 8, and 10 µm for exploring the G w -R w coupling effect.To reveal the effect of groove depth ( G D ), live cell experiments on G2NR, GNR2, and GNRN are performed and analyzed with three different groove depths (i.e., G D = 330, 725, and 1000 nm).

Cell morphology
In this section, three key cell morphological parameters (spread area ( A cell ), aspect ratio ( α ), and alignment angle ( θ )) of fibroblast cells are quantified by performing image analyses of live cells on a contact guidance (see Fig. 1a).
Our image analysis on the cell spread area indicates that cells on the patterned substrate tend to have smaller spread areas compared those on the flat control substrate (see Fig. S2 in the Supplementary Document).As a reminder, image analysis techniques are discussed in "Materials and methods" section.The spread area ( A Cell ) on the control substrate (= 2816.8 µm 2 ) is significantly larger than that on the topographical cues, but the variation in A Cell among different topographical cues is statistically insignificant ( A Cell = 2365.3,2115.3 and 2259.4 µm 2 on G D = 330, 725, and 1000 nm, respectively).Consistent with our results, a similar trend was observed for NH3T3 fibroblast cells cultured on micropatterned polydimethylsiloxane (PDMS) substrate 39 .
Figure 2 shows the elongation of cells.More specifically, Fig. 2a summarizes the cellular aspect ratio ( α ) as a function of N (the size of lateral widths) for three different cases: G2RN ( G w = 2 µm and R w = N µm ), GNRN ( G w = R w = N µm ), and GNR2 ( G w = N µm and R w = 2 µm ).In addition, different groove depths ( G D = 330, 725, and 1000 nm ) are considered.Note that Cont.indicates control experiments on flat control substrates.For clarification, the results of cell alignment and cell speed will be presented in the same way in the later sections to be consistent.
Consider the cellular aspect ratios from G2RN (see the first column in Fig. 2a).Two general trends are (1) α decreases with increasing R w and (2) α increases with increasing G D .For instance, α values of G D = 725 nm decrease from α = 13.06 (G2R2) to α = 5.81 (G2R10).Also, the distribution of the α values is sensitive to G D , e.g., a rate of decreasing α with increasing R w is much faster for G D = 725 and 1000 nm compared to G D = 330 nm.The strong G D dependency in G2RN can be again validated by comparing the average aspect ratio ( α ave ) for each G D .For example, α ave = 8.66 and 7.64 for G D = 725 and 1000 nm, respectively, are significantly larger than α ave = 5.92 for G D = 330 nm and α ave = 4.81 for the flat substrate (control).Kruskal Wallis test (see Fig. S3 in the Supplementary) confirms that the effect of G D on α is statistically significant (p << 0.001).The results of GNRN (second column) and GNR2 (third column) show somewhat similar qualitative trends of G2RN.However, our statistical test shown in Figs.S3 and S4 reveals that the p-values for varying R w are much smaller than those for G w regardless of G D .The smaller p-values are mainly caused by more statistically significant pairs in each case (e.g., G2RN and GNR2) rather than small variance within each topography by comparing standard deviations of cell aspect ratio for varying R w and G w .Therefore, we conclude that R w has stronger effect on the cell aspect ratio compared to G w .
It is important to acknowledge that α values depend on cell type as well as cell culture conditions.For example, different α values have been reported for cells on flat control substrates, e.g., α = 2-2.3for NIH 3T3 fibroblasts and myofibroblasts 39,40 .To account for a natural tendency of such cell-type and cell-culture dependent behavior, here all the aspect ratios are normalized by the aspect ratio of cells on the control substrate (i.e., α = α/ α Cont ).Following the previous studies [39][40][41] , the percentage population ( P α≥2 ) of elongated cells are quantified by using α ≥ 2.
Figure 2b shows P α≥2 values as a function of dimensional parameter N for cells on (left to right) G2RN, GNRN, and GNR2.In each plot, black diamonds, blue asterisks, and red triangles indicate G D =330, 725, and 1000 nm, respectively.First, P α≥2 for G D = 330 nm is much less sensitive to N compared to the deeper groove depths.As an example, due to the change in N from 2 to 8 µm , P α≥2 of G2RN decreases from 86.7 to 16.7% on G D = 725 nm and from 76.7 to 13.3% on G D = 1000 nm, respectively.The similar trend can be seen in GNR2 as the corresponding P α≥2 values are 86.7% ( G D = 725 nm) and 76.7% ( G D = 1000 nm) at N = 2 µm and 43.3% and 16.7% at N = 8 µm , respectivtly.Clearly, P α≥2 for G D = 725 and 1000 nm decreases with an increase of N at a much faster rate compared to G D = 330 nm, which indicates a strong depth dependent cell behavior.To further consider the correlation between cell elongation and groove depth for all lateral dimensions (R w and G w ), the average P α≥2 values for each depth, including all G2RN, GNRN, and GNR2 cases, are shown in see Fig. S5a.The average P α≥2 values are 4.44%, 16.5%, 31.9%, and 25.0% on G D = 0 (control), 330, 725 , and 1000 nm, respectively.We quantitatively characterize the effect of individual topographical features (R w , G w and G D ) and their products on P α≥2 by conducting multiple linear regression analyses.Because P α≥2 strongly depends on G D and changes its trend between G D = 725 nm and = 1000 nm, we performed the linear regression analyses twice using G D = 0−725 nm (Case 1) and G D = 0−1000 nm (Case 2), respectively, as summarized in Table S2 (see the Supplementary Materials).For clarity, we use G D , G w , and R w = 0 for the flat substrates (control) as no topo- graphical effect is expected on cell responses and, as a result, a constant term in the multiple linear regression model (see Eq. S1) becomes the P α≥2 value for the control substrates.A few key conclusions are (1) R 2 values for Case 1 and 2 are significantly different (i.e., 0.835 for Case 1 and 0.612 for Case 2, respectively) and (2) despite the different R 2 values, the coefficient associated with G D , for both cases, is 1-2 orders of magnitude larger than other coefficients with statistical significance (p < 10 −8 ).The smaller R 2 value for Case 2 may confirm a change in the P α≥2 trend for G D > 725 nm.
Finally, Fig. 3     www.nature.com/scientificreports/For more detailed analyses, Fig. 4a shows the cellular alignment angle (θ) for G2RN, GNRN, and GNR2.Note that θ increases with either the increase of G w and R w or the decrease of G D .For G D = 725 nm (the sec- ond row), the θ ave values for G2RN, GNRN, and GNR2, increase from 2.87 • ( N = 2 µm ) to 17.15 • , 11.74 • , and 4.69 • ( N = 8 µm ), respectively.Our statistical analysis in Figs.S6 and S7 indicates that the effect of G D is more significant than R w and G w in the light of the smallest p-value caused by more statistically significant pairs, not small variance within each topography.Also, it appears that the relative influence of G w on θ, compared to the other dimensional cues, is the weakest by the comparison of the corresponding p-values.
Figure 4b shows the percentage population of aligned cells ( P θ ≤2 • ) for different G w , R w , and G D values.It is worth noting that we have introduced a new criterion ( θ ≤ 2 • ) to quantify aligned cells because of the following reasons Other prior studies 27,42,43 categorized θ = 0-15° as aligned cells, but their standard is not adequate for our study since > 70 % cells on our patterned substrates are categorized as θ ≤ 15 • .On the contrary, we found that about 23% of cells on our chips satisfy θ ≤ 2 • , which also corresponds to the percentage population of elongated cells in Fig. 2b.Note that our new criterion allows us to consistently analyze the trends of P θ ≤2 • using the com- parable size in subpopulations as P α≥2 and, arguably, the correlation between P α≥2 and P θ ≤2 • can be considered.
For G D = 1000 nm (red triangle), P θ ≤2 • decreases from 50% ( G2R2 ) to 23.3% (both G2R8 and G8R2), whereas the corresponding P θ≤2 • values for G D = 330 nm (black diamond) are 13.3% ( G2R2 ), 6.67% ( G2R8 ), and 16.7% ( G8R2 ).G w and R w effects are more noticeable when G D values are 725 and 1000 nm.In addition, P θ ≤2 • is sensi- tive to G D up to G D = 725 nm (p < 0.001), but no statical difference between G D = 725 and 1000 nm (p ~ 0.583) as indicated in Fig. S5b.As an example, the average P θ ≤2 • values for G D = 0, 330, 725, and 1000 nm, regardless of G w and R w , are 3.33%, 10.6%, 29.8%, and 33.8%, respectively.It is worth noting that we performed regression analyses of P θ≤2 • using Eq.S2 and obtained similar conclusions (see Table S3 in the supplementary) compared to the regression results of the P α≥2 data.The effect of groove depth on the cell alignment is the most significant cue.
α and θ are the most sensitive to G D among the controlled morphological cues in this study.To further con- sider the effect of G D on both α and θ , regardless of G w and R w , Fig. 5 correlates the cell alignment angles with the corresponding cell aspect ratios for different groove depth ( G D ).Interestingly, the α-θ correlation appears to be inversely proportional.For example, for G D = 0(control), 330, and 725 nm , the trend is a decrease in θ with an increase in α .Then the direction of the trend changes from 725 to 1000 nm as evidenced by a similar θ (p ~ 0.44), but a smaller α (***p < 0.001) with the considerable statistical significance.
As mentioned above, there is a significant discrepancy in the previous reported and our new θ values likely due to cell-type-dependent contact guidance.Moreover, Crouch et al. 27 studied human foreskin fibroblasts (HFF) only using G5R5 configuration with different depths from G D = 800 nm to G D = 1024 nm .Despite the nearly identical cell aspect ratio of HFF and Hs27 cells (Fig. S9) as a function of the groove depth, the corresponding cell alignment again exhibits considerable cell type dependent discrepancy.Because of these cell type dependency, time-effective characterization of cellular response to contact guidance, e.g., using multiplexed contact guidance chips to concurrently consider a wide range of topographic design parameters within the same cell culture, becomes increasingly relevant for practical optimization and application of contact guidance technologies for biomedical applications and biological studies.

Cell migration: dynamic image analysis
In this section, we perform analyses to quantify dynamic cell migration on contact guidance topographies with an emphasis on directional migration along the groove direction.All the directional characteristics of the dynamic cell migration in the following analyses are characterized with respect to the groove direction.
Figure 6a and Fig. S10 show the overall migration directionality.For this, vectors from the initial to final positions of the same individual cells migrational tracks > 10-h intervals were obtained from optical images.Then the angle between the vector and the groove direction was calculated for the migration directionality.Note that the percentage cell population near 0 • monotonically increases with the groove depth.For example, the cell populations between 0 • and 10 • are 8.21%, 41.7%, 71.8%, and 75.3% on G D = 0, 330, 725, and 1000 nm, respec- tively.For statistically significant conclusions on the G D -sensitive directional orientation, we have performed Kruskal Wallis test as shown in Fig. S11a.The average directional orientations are 43.8 • , 21.6 • , 10.7 • , and 9.83 • on G D = 0, 330, 725, and 1000 nm, respectively.The results confirm that the directional orientation decreases from G D = 0 to = 725 nm ( p ≪ 0.001 ) and then the difference between G D = 725 and = 1000 nm becomes statistically insignificant (p ~ 0.06).It is worth noting that the effect of R w on the cell directional orientation is statistically obvious (p < 0.05) only for G D = 330 nm among three different depths (see Fig. S11b,c).Moreover, for G D = 330 nm, R w has a stronger influence on the cell directional orientation than G w (i.e., R w : p ~ 0.0002 and G w : p ~ 0.063), which results from more statistically significant pairs instead of small variance within each topography group as shown in Figs.S11b,c and S12.
The cell directional orientation provides the general direction of cell migration over > 10 h but misses timedependent dynamic characteristics in a shorter period.Because of this limitation, we compute an angular displacement (dθ) of individual cells using different time intervals (i.e., δ = 10, 50, 100 and 200 min) and map the occurrence distribution (i.e., the normalized cell population) with a 15 • -bin size using the following expression 44 : where − → v t is a velocity vector of a cell between t and t + δ over a finite time interval ( δ ) (see the inset in Fig. 6b).The occurrence in Fig. 6b exhibits two general trends for all three G D cases.First, the occurrence for δ = 10 min exponentially decays with increasing dθ , but a larger G D value results in faster decay.Second, the occur- rence in 0 • ≤ dθ ≤ 15 • decreases, while the occurrence in 165 • ≤ dθ ≤ 180 • increases, with increasing a time interval, respectively.These observed trends indicate that cell migration on deeper G D is indeed more direc- tional in the groove direction for any time interval.For example, note that the occurrence in 0 • ≤ dθ ≤ 15 • on G D = 725 and = 1000 nm (e.g., 0.285 and 0.303, respectively, for δ = 200 min) is significantly higher compared to G D = 0 nm (e.g., 0.276 even for δ = 10 min).In addition, the occurrence in 165 • ≤ dθ ≤ 180 • for deeper G D increases with increasing δ at a faster rate compared to the control, for example, from 0.04 ( δ = 10 min ) to 0.16 ( δ = 200 min ) for G D = 1000 nm while the occurrence for the control remains low (i.e., 0.09 for δ = 200 min ).In summary, cells on deeper G D exhibit strong forward motion at shorter δ and back-and-forth motion at larger δ , respectively, along the groove patterns.Our additional analysis, while not discussed in detail here, indicates that the influence of R w and G w on dθ is not significant (see Fig. S13a,b, respectively).
We next identify how individual topographical cues affect migration speed of Hs27 cells.Surprisingly, the cell migration speed difference between G2RN, GNRN, and GNR2 is insignificant and the average migration speed ( v = total travel distance over total travel time) is insensitive to R w and G w .Again, v strongly depends on G D (see Fig. S14).Because of this reason, we summarize cell migration speed for different groove depth, G D , in Fig. S15 without separately considering R w and G w .The mean value of the cell speed on the flat substrate (control) is about 0.4 µm/min , which matches with the previous study using human skin fibroblasts 45 .The migration speed on the G D = 725 nm is 0.25 µm/min , which is significantly lower than 0.39 and 0.35 µm/min on the control (1) www.nature.com/scientificreports/and G D = 330 nm , respectively.The average cell speed on G D = 1000 nm is 0.29 µm/min , slightly higher than G D = 725 nm.Interestingly, the observed trend in the average cell speed is qualitatively similar to the cell align- ment angle in Fig. 5.
To quantify the directional speed, the average migration speed is decomposed into parallel ( v y ) and perpen- dicular ( v x ) directions with respect to the groove direction as shown in Fig. 7a.The v x and v y are calculated by total travel distance along x and y directions over total travel time, respectively.The v y to v x ratio increases from 1 (control) to 3.1 ( G D = 725 nm ) with an increase of G D , but it slightly decreases to 3.0 with further increases of G D = 1000 nm .Note that v y is about 3.0-3.1 times larger than v x for G D = 725 and 1000 nm , indicating direc- tional cell migration in the groove direction.In contrary, v x ≈ v y for the control experiment indicating that the cells move without a preferred direction.
So far, we have quantified the average cell speed under the influence of contact guidance.However, it is still crucial to amplify the groove direction motion (along y-axis) of cells by minimizing their lateral motion (width direction along x-axis) to effectively guide cells to target locations, e.g., injured areas, using the contact guidance technology.For an increase of migration efficiency, we consider the concept of y-axis directional persistence through a directional displacement ratio.In this section, our emphasis is placed on the effect of G D because the effect of G w and R w on cell speed is insignificant as discussed above (also see Fig. S16).
The question to be addressed is whether as an efficient directional guidance, the cell migration trajectory can be effectively controlled by a change of G D scales in order to rapidly reach our target point (i.e., somewhere along the groove direction).Similar with well-known directionality ratio (i.e., d D , d = straight-line length between the start point and the endpoint of the migration trajectory, and D = the total length of the trajectory) 46 , the directional displacement ratio has been utilized to characterize y-axis directional efficiency, which is defined as follows: In other words, instead of using d, the directional displacement is the ratio of a projection of d onto the y axis ( d y (0, �) ) over the total length of the trajectory ( D(0, �) ) from start ( t (time) = 0) to finish ( t = ) in which the time interval ( � (min) ) ranges from 10 to 600 min with 10 min increment as shown in the inset of Fig. 7b.Since the d y and D depend on the , this ratio can markedly change over the time period of a cell trajectory.The directional displacement ratio for all cases in Fig. 7b decays exponentially in , but a larger G D value result in slower decay.More specifically, the initial values of the directional displacement ratio are 0.65 (control), 0.72 ( G D = 330 nm ), 0.78 ( G D = 725 nm ), and 0.78 ( G D = 1000 nm ) and then they decrease to 0.26, 0.43, 0.54, and 0.48, respectively, at = 600 mins.It is worth remarking that the directional displacement ratio for the control is consistently smaller than those for all three G D cases as expected.The important observation, here, confirms that the cells require less total length of the trajectory to arrive at a target location along the groove direction with a deeper G D .

Discussion
Quantitative characterization of cell morphologies associated with a wide range of microfabricated surface patterns is relevant to understanding and guiding cell behavior for several biological applications.For example, cell spread areas have been studied in the context of cell differentiation and proliferation 47 , migration 48,49 , and transfection 50 .Note that the cell spread area becomes important for fibroblasts as it regulates cell migration in wound healing process 51 .It is also known that cell shape, often quantified by a cellular aspect ratio, is related to regulation of division-coupled interspersion for cells 39,52 .Lastly, aligned organization of cells, also known as a cell polarity, contributes to tissue maturation 53 , regeneration growth 54 , differentiation 55 , and proliferation 56 and, as a result, a capability to control cell alignment is crucial for engineering functional tissue equivalents.Accurate characterization of dynamic cell behavior on various topographical cues is crucial for many biological applications such as embryological development, tissue formation/regeneration, immune defense/inflammation, and cancer progression [57][58][59] .For instance, it is known that the directional cell movement of fibroblast enables development of a larger contact area during wound healing process 60 .Furthermore, cell migration speed is directly linked to acceleration/deceleration of pathological matrix deposition processes in fibrosis 61,62 .Finally, an overall efficiency of cell migration [63][64][65][66][67] , typically quantified by the directional motion and its ratio over the total travel distance, is relevant to rapid ECM synthesis and remodeling.Our study may offer a means to efficiently guide cell migration in the scope of achieving a quicker recovery of the microenvironments 68 .
In this study, we have characterized the influence of different topographical cues (i.e., R w , G w , and G D ) on the static and dynamic behavior of Hs27 fibroblast cells by utilizing a microfabricated contact guidance chip.To map the effect of the R w and G w beyond what has been reported in the literature, we have considered three different cases, namely, G2RN, GNRN, and GNR2.Our study on the cell morphologies shows that the cells are likely elongated and aligned with the groove direction as G D increases from 0 to 725 nm with the correspond- ing p-value < 10 -8 .Within the same G D range, the effect of R w is also statistically significant (p-value < 0.005).It appears that the trend in cell responses changes between G D = 725 and 1000 nm because the cell elongation is statistically similar while its alignment slightly decreases as evidenced by a considerable decrease in the R 2 value.
Our experiments conclude that the groove depth, G D , is the most dominant topographic cue for both static and dynamic responses of Hs27 fibroblasts, which agrees with other studies 11,27,69 .To provide feasible explanations of the G D -dependent cell behaviors, we discuss the role of cell membrane on contact guidance.When G D increases for a given ridge and groove width, cell membrane must undergo a larger degree of membrane deformation to occupy the space within the cross section of individual grooves defined by G w and G D .A cell membrane protrusion led by filopodia will sense topographic cues and move toward the groove bottom.Considering that the average length of filopodia is about 1-5 μm [70][71][72] , the Hs27 cells will likely reach the groove bottom and form focal adhesions as confirmed in our previous study 35 .For simplicity, we assume that the lateral filopodia protrusions occur on the ridge edge with their characteristic length ( l f ).Then, the bending angle ( θ f ) between filopodia and the groove bottom can be calculated as follows: The Eq. ( 3) shows that θ f increases with increasing G D .In Ref. 73 , it was shown that larger θ f increases the normal traction force component on the integrin-based filopodia adhesion.Such high normal stresses lead to an impairment of FA formation as well as severe buckling of the branched actin fibers in the filopodia 70,73,74 .As indicated by several studies [75][76][77][78][79][80] , buckled F-actin filaments shorten and break the actin shaft and, therefore, filopodia forms in the lateral direction will likely be short-lived and retracting.Moreover, Dunn et al. reported that a larger bending angle, e.g., to follow discontinuous topographic patterns on a contact guidance chip, further limits lateral protrusions of filopodia across the features 3 .Based on the previous studies, filopodia likely becomes thick, actin-rich, and stable in the groove direction compared to the lateral direction 81 .This may offer possible explanation for the directional cell orientation and migration because filopodia is known to be linked with the cell morphology 82 as well as the cell migration 83,84 .
For the R w -dependent morphology of Hs27 fibroblasts, a model based on the constrained maturation of focal adhesions (FAs) 16,85 can provide mechanistic insights.The anisotropic and oriented topographical cues, i.e., parallel ridges and grooves, likely induce cell membrane protrusions as well as FA formation and maturation along the topography.These laterally restricted FAs result in the formation of aligned actin stress fibers to along the parallel direction of the topography because the orientation of the stress fibers follows that of the FAs 35 .Note that the aligned actin stress fibers connected to large, long-lived, and oriented FAs are longer and more stable than small, short-lived, and non-aligned ones associated with non-oriented FAs 16 .With reduction of R w , the formation of FAs is more confined along the groove direction and, as a result, the Hs27 fibroblasts tend to be more elongated and aligned in the same direction as shown in Fig. 4.
Our dynamic analyses also confirm G D as the main geometrical parameter to influence cell behavior responses.The cells on deeper G D are better guided along the groove direction.Cells on the flat control and G D = 330 nm substrates move faster, compared to those on G D = 725 and 1000 nm .This trend ( v = 0.55 µm/min for control and v = 0.36 µm/min for the grooved pattern) has also been observed by other study which uses NIH 3T3 fibroblasts with R w = G w = 0.43 µm and G D = 200 nm 86 .However, these results are directly opposite with those of another study which reports that faster migration speed for MC3T3-E1 mouse osteoblast cells on various G D ranging from 0 ( v = 0.39 µm/min ) to 4.5 µm ( v = 0.61 µm/min ) with R w = 6 µm and G w = 4 µm 20 .Based on a positive feedback loop between actin flow rates and the cell elongation, the elongated cells are likely to accelerate their actin retrograde flow such that it improves their motility as confirmed using several different cell types such as retinal pigment epithelial cells and bone marrow-derived mice dendritic cells 15 .This discrepancy may come from not only cell type-, but substrate topography dependent characteristics.

Conclusion
With growing interests in the development of the next generation biomaterials, a fundamental understanding of contact guidance offers the potential to guide specific cellular functions in vivo.Our current work characterizes static and dynamic behaviors of the Hs27 fibroblasts in the wide range of groove and ridge dimensions and offers better understanding of individual geometric factor's contribution on the cellular responses.For (3) static results, groove depth ( G D ) as deep as 725 nm induces significant elongation ( α ave = 8.66 ) and alignment ( θ ave = 7.23 • ), while further increase of G D to 1000 nm slightly decreases both elongation ( α ave = 7.64 ) and alignment ( θ ave = 9.27 • ) compared to the cells on a flat substrate (control, α ave = 4.81 and θ ave = 39.7 • ).The regression analyses show that G D and G D • R w terms have a statistically significant effect on the cell elongation and alignment, whereas the effect of G D as deep as 1000 nm has been diminished compared to that of G D = 725 nm based on its coefficient reduction.In our dynamic analysis, a deeper G D better guided the cells migrating along the groove direction, even though their average speed slightly decreases from 0.39 µm/min (control) to 0.29 µm/min (1000 nm), respectively.The high occurrence of the 0 • ≤ dθ( angular displacement) ≤ 15 • and 165 • ≤ dθ ≤ 180 • from the groove direction and its faster decaying rate with increasing dθ indicate the y-axis directional movement of the cells on the deeper G D .The large ratio of directional speed and -directional displacement, both in the groove direction, for a deeper G D suggest that the cells on the deeper G D are better guided by minimizing their unwanted motion perpendicular to the groove direction.These findings demonstrate the importance of G D for the Hs27 fibroblasts as a guideline of the efficient migration journey.Our in-depth experimental and statistical studies would provide a way of creating an effective platform for the application of topographical cues, which are shown to be successful in guiding and controlling cell morphological and migrating characteristics.
Our multiplexed, unidirectional groove structures enable us to systematically investigate the effect of each physical constraint on the cellular responses with tightly controlled other substrate features including surface chemistry and mechanical properties (e.g., stiffness).The microfabricated topographical features aim at mimicking physiological topography such as aligned fibers and ECM tracks that cells frequently encounter in vivo.With the use of our topography, we can design well-controlled and simplified in vitro systems that allow us to tightly control cell static and dynamic behaviors as close as those on the natural surroundings.This attractive feature may offer unique research opportunities to simulate cell morphology dependent behaviors, e.g., how cell membrane damage is linked to cell morphology.

Figure 1 .
Figure 1.(a) A monolithic quartz contact guidance platform (left) and a schematic diagram of how cell topography is characterized on the micropatterned substrate (right).A-A is a cross sectional view of the platform showing groove (G w ) and ridge (R w ) widths, and groove depth (G D ).The outline of a cell is fitted by an ellipse to estimate a length ratio of the major to minor axes and a cell alignment angle ( θ ) between the major axis and the direction of the groove patterns.(b) Optical images of live Hs27 cells on G w = R w = 10 µ m (i.e., G10R10) with (left) G D = 330 nm, (middle) = 725 nm, and (right) = 1000 nm, respectively.
13:21691 | https://doi.org/10.1038/s41598-023-48913-9 shows the normalized cell population as a function of the cell alignment angle (θ) on the control and patterned substrates with three different groove depths ( G D ).Despite significantly different population dis- tributions shown in the figure, the average values of the cell alignment angle ( θ ave ) for G D = 0(control), 330, 725, and 1000 nm are 39.68 • , 19.70 • , 7.23 • , and 9.27• , respectively.For clarification, each average value for each G D includes G2RN, GNRN, and GNR2 and, therefore, the effect of G D is considered for all R w and G w .Note θ for the patterned substrates are significantly smaller than that of the control substrate and their distributions are also statistically different with p-value << 0.001 as shown in see Fig.S6.The cells on the control substrates exhibit a uniform distribution of alignment angle distribution unlike the patterned substrates.This expected trend can be also seen in the linear increasing cumulative distribution function (CDF) plots shown in Fig.S8.

Figure 3 .
Figure 3. Normalized cell population as a function of cell alignment degree with respect to surface patterns: (from left to right) control, 330, 725, and 1000 nm.The bin size is 2 • .

Figure 5 .
Figure 5.The cell alignment ( θ ) as a function of the cell aspect ratio ( α ) for different groove depths.The error bar denotes 95% confidence level interval.

Figure 6 .
Figure 6.(a) Directional orientations of the fibroblasts on contact guidance with different groove depths.The angular and radial coordinates represent the cell orientation from 0 to 90 • and the corresponding percentage cell population within the total population (each bin size is 2 • ).(b) The percentage population distribution of angular displacement ( dθ , see graphical definition in inset and δ represents a time interval) evaluated at four different time intervals (10, 50, 100, and 200 min).For (a,b), each column represents the results of different groove depths (left to right) from 0 to 1000 nm.

( 2 ) 2 .Figure 7 .
Figure 7. (a) Cell migration speed in the groove direction ( v y ) and perpendicular to the groove direction ( v x ).(b) Directional displacement ratio (i.e., y direction moving length ( d y (0, �) ) divided by the total travel length ( D(0, �) ) from start (t (time) = 0) to finish ( t = ), see graphical definition in inset) as a function of the time interval ( ) for different groove depths.All error bars are mean ± standard error of mean (SEM).